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Abstract — In this paper we develop a low-complexity coding 
scheme and system design framework for the half duplex relay 
channel based on the Quantize-Map-and-Forward (QMF) relay- 
ing scheme. The proposed framework allows linear complexity 
operations at all network terminals. We propose the use of binary 
LDPC codes for encoding at the source and LDGM codes for 
mapping at the relay. We express joint decoding at the destination 
as a belief propagation algorithm over a factor graph. This graph 
has the LDPC and LDGM codes as subgraphs connected via 
probabilistic constraints that model the QMF relay operations. 
We show that this coding framework extends naturally to the high 
SNR regime using bit interleaved coded modulation (BICM). We 
develop density evolution analysis tools for this factor graph and 
demonstrate the design of practical codes for the half-duplex 
relay channel that perform within IdB of information theoretic 
QMF threshold. 

L Introduction 

Cooperative relaying has been proposed as a promising 
technique to resolve the increasing demand for data throughput 
in wireless networks. Recently a lot of progress has been 
made in establishing the theoretical foundations of cooperative 
communication. To apply these principles towards the design 
of practical wireless systems, various system design tradeoffs 
must be taken into consideration. This paper presents progress 
towards this goal. We propose a system design and coding 
framework for quantize-map-and-forward (QMF) [IJ relaying 
that has low complexity and performs close to information 
theoretic bounds. 

A. Cooperative Systems 

A cooperative wireless link typically consists of an infor- 
mation source, a destination and one or more cooperating half 
duplex relays. The relays are usually assumed to operate in- 
band. i.e. no additional channel resources are allocated for 
cooperation. Without loss of generality, it is assumed that re- 
lays use time-division-duplexing i.e. they listen to transmission 
from the source for some fraction of total time, then forward 
a description of their observation in the remaining fraction. 

There are several aspects involved in the design of a 
cooperative relaying system. Listening fractions and forward- 
ing schemes must be determined for each relay. Suitable 
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modulation and channel coding schemes must be designed 
for various terminals. Rate adaptation mechanisms must be 
considered to account for changes in availability of relays and 
channel strengths. Practical constraints must be considered e.g. 
minimizing the overall system complexity, reuse of building 
blocks from traditional (non-cooperative) systems as much 
as possible, compatibility with protocols at higher layers 
and handling of system imperfections like synchronization, 
channel estimation errors etc. In this paper, we focus on the 
coding and signaling aspects of cooperative relaying. Other 
components are discussed in brief towards the end of the paper 

B. Relaying Schemes 




Fig. 1. Example relay network. With multiple antennas at destination, source- 
relay cooperation provides additional degrees of freedom for communication. 

Most wireless systems operate at moderate to high SNR i.e. 
in a regime where transmit power is not the major limiting 
factor for link capacity. At high SNR, the link capacity is 
limited by spatial degrees of freedom. Relay cooperation is 
of special interest for practical systems because it has the 
potential to provide additional spatial degrees of freedom. 
For illustration, consider a relay channel with single- antenna 
source, a half-duplex single antenna relay and a destination 
with two antennas shown in Fig. [T] Since the destination has 
multiple antennas, the source can spatially multiplex traffic to 
the destination using relay cooperation. If the source-to-relay 
channel is strong, this network can approach the high- SNR 
performance of the 2 x 2 MIMO channel O. 

Various strategies for relay cooperation are proposed in 
literature. Among these amplify-and-forward (AF), decode- 
and-forward (DF) and compress-and-forward (CF) O, lU 
have received the most attention. Under DF, the relay decodes 
the source's message and forwards a hard estimate of it, 
whereas under AF and CF it forwards a soft estimate without 
explicitly decoding it. In DF and CF, the relay maps its 
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estimate to a random codeword before forwarding, whereas 
in AF it forwards an uncoded signal. The QMF scheme |T| 
also uses soft estimate forwarding with random coding similar 
to CF. For the example network in Fig. [T] the CF and QMF 
schemes are close to optimal at high SNR. In fact they achieve 
within one bit/sec/Hz to the information-theoretic capacity |5 |, 
im . An intuitive explanation for why QMF/CF performs better 
than both AF and DF is given in the context of the example 
in Fig. [T] below. 

In the example network of Fig. [T] the destination receives 
continuously from the source. The relay receives a stronger 
version, when it is listening. Since the destination has two 
antennas, it can resolve simultaneous transmissions from the 
source and relay. In order to achieve spatial multiplexing, 
the relay should extract the less significant bits (from its 
observation), which the destination cannot resolve, and for- 
ward them. Under AF, the relay only forwards the more 
significant bits, which the destination can already resolve. 
Therefore AF cooperation provides limited benefit. Under 
DF, the relay decodes the entire message before it forwards 
anything. Since the listening time is limited, this approach is 
inefficient. The QMF/CF schemes implicitly extract the less 
significant bits from the relay's observation by using quan- 
tization/compression and random mapping. Therefore these 
schemes provide the most cooperation gain. 

Despite having similar performance for the single relay 
network, the CF and QMF schemes have significant differ- 
ences. In the conventional CF scheme, the relay compresses 
its observed signal and performs a random code mapping 
before forwarding. The compression rate is chosen in order 
for the destination to perform two-step decoding i.e. first 
decode compressed signals from relay and then use it as 
side information to decode the message from source. For 
configurations that involve multiple relays, two-step decoding 
is sub-optimal and conventional CF is not within bounded gap 
from information-theoretic capacity [IJ. Even for the single- 
relay configuration, conventional CF requires that the relay 
have full knowledge of the quality of its forward channel. 
This introduces a large estimation and feedback overhead for 
fading channels and increases the complexity of rate adaptation 
schemes. 

Under QMF, the relay quantizes its received signal at noise 
level, randomly maps it to a codeword and forwards it. Unlike 
CF, the quantization and mapping is performed without regard 
to the quality of forward channel at the relay. This reduces 
the channel estimation and feedback overhead for the link. It 
also simplifies rate adaptation protocols. Additionally, QMF 
uses joint decoding (as opposed to successive decoding) and 
performs within bounded gap from capacity for networks 
having an arbitrary number of relays [IJ. QMF has played 
a key role in several recent information theoretic results on 
cooperative networks 0, 01, ||2l, ||7J. Due to these favorable 
properties, the QMF scheme is superior to CF from the 
perspective of practical cooperative systems. 

Since mapping at a QMF relay is performed without any 
knowledge of forward channel strength, side information from 
relays cannot be decoded at the destination independently. 
QMF requires joint decoding of the message (from source) 



and side information (from relays) IH. This presents a unique 
challenge because joint decoding typically requires higher 
complexity and makes it harder to design a practical coop- 
erative coding scheme. The key contribution of this paper is 
to develop a low-complexity cooperative coding framework 
for QMF that significantly reduces the complexity of joint 
decoding and yet performs close to information theoretic 
bounds. 

C. Related Work 

Majority of previous work on code design for cooperative 
relaying is focused on the DF scheme. DF relays fully decode 
the source's message. Therefore, DF coding schemes involve 
partitioning a large codebook into two parts. The source 
transmits one part of the codeword and the relay transmits 
the remaining part IH, ||9l, MPI. Turbo code designs which 
perform ^ IdB away from the DF information theoretic 
threshold are demonstrated in lITTIl . lll2l . LDPC profiles are 
developed for DF in (TSl. A bilayer LDPC structure dl 
and the protograph method ifTSl has been used to get LDPC 
designs < 0.5dB from the DF threshold. The bilayer structure 
is extended for use at high SNR using bit-interleaved coded 
modulation (BICM) |[T6l . As for CF relaying, a coding scheme 
using a combination of LDPC and irregular repeat accumulate 
(IRA) codes is presented in 1 17 |. Rateless coding schemes are 
developed in ifTSll . As for QMF relaying, a coding scheme is 
proposed in independent work 1 19] based on lattice strategies. 
The scheme in |19 | reduces the complexity of mapping at the 
relay to polynomial-time while the joint decoding complexity 
remains exponential- time. 

D. Summary of Results 

In this paper, a coding scheme for QMF relaying with 
linear complexity encoding at the source, mapping at the 
relays and joint decoding at the destination is developed. For a 
network with one relay, the proposed scheme performs within 
(0.5— l)dB gap from the information-theoretic QMF threshold. 
For the code design example considered in Sec. |Vl the QMF 
threshold is ^ 1.5dB better than DF. The key techniques used 
in this paper are summarized as follows: 

1) BICM: Design of binary channel codes with standard 
higher order signal constellations is considered based 
on the widely used BICM technique O. 

2) LDPC-LDGM: The scheme uses low density parity 
check (LDPC) codes at the source for channel coding 
and low density generator matrix (LDGM) codes at the 
relays for mapping. 

3) Joint Factor Graph: The joint decoding procedure at 
the destination is formulated as a belief propagation 
algorithm over a factor graph. This graph contains the 
original channel code (LDPC) and relay mapping func- 
tions (LDGM) as subgraphs connected via probabilistic 
constraints that model the QMF relay operations. 

4) Practical Decoding Algorithm: Using a DBLAST space- 
time architecture, scalar quantization procedure at relays 
and specific choice of component codes, the resulting 
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factor graph is greatly simplified, making it suitable for 
practical decoder implementation. 
5) Code Design: Density evolution analysis tools 1211, (221 
are developed for the systematic design of joint LDPC- 
LDGM factor graphs. 

E. Organization 

In Sec. nil the coding framework for QMF and corre- 
sponding joint decoding algorithm is developed. The treatment 
focusses on a canonical system model with one relay and 
binary inputs. In Sec. [Till density evolution and code design 
tools are developed. In Sec. |lVl the framework is extended 
to the high SNR regime i.e. for high order modulation inputs 
using BICM. In Sec. |Vl the design of codes for an example 
cooperative link is demonstrated. Finally in Sec. |Vl] a sketch 
is provided for extending the proposed framework to scenarios 
with multiple relays. 

II. Coding Framework 

A. System Model 

Initially, this paper focuses on the design of codes for a 
binary memoryless symmetric (BMS) relay channel as de- 
scribed below. In Sec |lVl this model is extended to high order 
modulation inputs for use at high SNR. 

The BMS Gaussian relay channel has three half-duplex 
terminals: source {S), relay {R) and destination {D) with 
binary input additive white Gaussian (BIAWGN) channels 
between them, as shown in Fig. [21 R listens for a fraction 
/ G [0, 1] of the total communication time and transmits for 
the fraction (1 — /). The block lengths for the transmitted 
codewords at S and R are Ns and Nr respectively. They 
satisfy the half-duplex constraint Nr = (1 — f)Ns. The 
codeword messages sent from S and Rwcehs G {0, 1}^^ and 
bi? G {0, 1}^^ respectively. The corresponding transmitted 
signals are G {±^/Ps}^^ and x^^ G {±\^Pr}^^ where 
Ps and Pr are per-node symbol constraints on average power 
i.e. E\x\ ■\ < Ps and E\x\-\ < Pr. Bold-face lower case 
letters are used to denote a sequence of symbols. 




Fig. 2. Half-Duplex Binary Input Gaussian Relay Channel 

Multiple (M) receive antennas are assumed at the destina- 
tion. This permits consideration of network scenarios where 
cooperative spatial multiplexing is possible i23ll ll2ll. An ex- 
ample scenario with M = 2 is discussed in the Appendix. 
The received signals at D and R are denoted as G 
and yR^j for each symbol time i G {1,2,..., Ns} and j G 
{1,2,..., fNs} respectively. They are modeled as follows: 



Here hi,h2,/ii? denote the corresponding channel gains. 
x'^. = and h2 = for i G {1,2,... J Ns} when R 
is listening. For the remaining time x'^^- = XR^i-fNs^ 
i G {fNs + 1, . . . , Ns}. Zi and zrj are i.i.d. zero-mean 
Gaussian noise vectors with identity covariance matrices. All 
the channel observations at D are denoted by Y G C^^^^ 
i.e. Y = [yi y2 ... JNs]- Observations at R are denoted as 
jR G C^^^-^^^ i.e. jR = [yR^i ... Vrjns]- The channel 
is characterized by the following parameters: SNR^^^ = 
Ps\hR\\ SmsD = Ps\\hi\\^ and SNRrd = Pi?||h2|p. 

B. Quantize-Map-Forward Scheme 

The quantize-map-and-forward scheme [IJ is summarized as 
follows. S has a sequence of messages rrik G {1, . . . , 2^^^}, 
/c = 1, 2, ... to be transmitted. At both S and R, codebooks 
Cs and Cr are created respectively. S maps each message 
to one of its codewords and transmits it using Ns symbols 
resulting in an overall transmission rate of IZ. Relay listens 
to the first fNs time symbols of each block. It quantizes its 
observation at noise level i.e. the quantization distortion is 
equal to the noise power at the relay. Relay maps the quantized 
bits to a codeword in Cr. It transmits this codeword using 
(1 — f)Ns symbols. The destination D attempts to decode 
the message sent by S from received signals (Y). In order to 
decode, D must know all channel parameters SNR^i), SNRi^i^) 
and SNR^i^, the relay listening fraction / and both codebooks 
Cs and Cr. 

It is assumed that SNR5'i),SNRi^i) are measured at D and 
SNR^i^ is measured at R using pilot symbols. It is further 
assumed that SNR^i^ is forwarded to D by R. The estimation 
and forwarding overhead of these steps is ignored for the 
analysis presented in this paper. 

C. Factor Graph for Joint Decoding 

In the context of the system model and cooperation scheme 
outlined above, let us focus on binary linear codebooks C| and 
Cj^. These can be represented as bipartite Tanner graphs using 
respective parity check matrices. In such a representation bit 
(variable) nodes represent the codeword and check (function) 
nodes represent parity constraints that must be satisfied in 
order for the codeword to be valid. Let us consider the 
maximum a posteriori (MAP) rule for joint decoding at D. 
In this subsection, joint decoding is expressed as a sum- 
product algorithm over a factor graph that contains the Tanner 
graphs of component codes (C^,Cj^) as sub-graphs connected 
via probabilistic constraints that represent the QMF relaying 
operation [24J|25||26|. 

Joint decoding involves searching for the codeword G 
Cg that maximizes the a posteriori probability p(b5'|Y). An 
efficient way to do this search is to consider the bitwise maxi- 
mum a posteriori (MAP) decoder, where the aim is to compute 
P{bsA^) = E^6.,.P(b5|Y) foralH = l,2,...Ar^. 

f{Y\hs,hR)p{hs,hR) 



p{hs\Y) = ^- 



/(Y) 



^ixs,i 



Vrj = hRXsj + zr^. 



^^f{Y\hs,hR)p{hs,hR). 
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For the first fNs bits, R is Hstening and D observes an 
interference-free signal from S. During the remaining trans- 
missions, D observes a superposition of signals from S and R. 
Therefore, the first term / (Y|b5',bi^) factorizes as follows: 

f{Y\hs,hn) 

fNs Nr 

i=l j=l 

The codes C| and Cj^ have characteristic functions l(b5' G 
C|) and l(bi^ G C^) respectively. 

p {hs, hn) = p (hs) p (b^^lb^) 

^l{hseC's)p{hR\hs) 

1 {hs e C's) 1 {hR e C'j,) p{hn\hs) . 
(a) is due to the fact that b^^ must be a codeword in Cj^. 

CHK P(bi?|bs) CHK 




/(yi|^5,i) fiy2\bs,2) \ X/ / 

OBS ■ ■ 

f{yifNs + l)\bs,ifNs + l)^bR,l) 

Fig. 3. Factor Graph for Joint Decoding. 



f{yNs\bs,Ns,bs,Nn) 



The resulting factor graph in Fig. [3] shows that in addition 
to nodes representing channel observations i.e. /(yi|6s',i, bRj) 
the subgraphs 1 (b^ G C|) and 1 (b^^ G Cr) are connected by 
p{hR\hs) that represents the quantization operation at R. 

If the component codes C| and Cr are sparse, the overall 
factor graph is also sparse. A sum-product algorithm for 
decoding over such a factor graph has complexity that grows 
linearly with the length of component codes. However, the 
sum-product update rules at the function node ^(bi^lb^) is 
very complex due to its high degree (Ns + Nr). Moreover, it 
introduces very short cycles in the graph, which deteriorates 
the performance of sum-product decoding. In order to get 
reasonably close to MAP performance and low decoding 
complexity, the p{hR\hs) node must be factorized further. 
In the following subsections, choice for component codes C| 
and Cr and specific techniques for factorization are discussed. 

D. Choice of component codes 

In the discussion above, general binary linear codes and 
C\ are considered. A natural choice is to use sparse graph 
codes (like LDPC) that are known to have good performance 
and linear complexity decoding/encoding operations. 

In a previous communication |24|, preliminary results for 
such factor graphs were presented using off-the-shelf LDPC 
codes at both S and R. As observed in ||24| . off-the-shelf 



(point-to-point) LDPC codes do not allow close-to-optimal 
performance over cooperative channels. An information- 
theoretic understanding of this observation is presented in 
|27|. The authors point out that capacity-achieving codes for 
the point-to-point channel exhibit higher estimation errors 
whenever the SNR is below the Shannon limit. Therefore, in 
cooperative networks where the operating SNR is below the 
point-to-point Shannon limit, such off-the-shelf codes are no 
longer suitable to utilize side information from the relay at the 
destination. As a consequence, specialized codes are required 
for cooperative channels. For sparse graph codes, specialized 
code profiles that are optimized for relaying can be designed 
using standard tools such as density evolution analysis EH, 
L22J. However, for the LDPC-LDPC combination fSt] density 
evolution does not extend readily to QMF joint factor graphs. 

In this paper, the use of LDPC codes at 5* and LDGM codes 
at R is proposed. LDPC codes are known to perform very 
close to information theoretic limits when used for channel 
coding. Similarly LDGM codes are commonly used for lossy 
data compression f28l and the LDPC-LDGM combination 
is a good fit for the QMF relay channel. Moreover, density 
evolution analysis tools can be extended to LDPC-LDGM joint 
factor graphs. Such an extension is developed in Sec. |llll This 
permits explicit construction of code profiles optimized for 
relaying. 

Based on the LDPC-LDGM choice, let us introduce aux- 
iliary variable nodes bg = {hQ^i}f^^ in the factor graph, 
bg represents the Kr bits after quantization at R. These 
are mapped to the codeword bR of length Nr obtained after 
passing through a low density generator matrix having Kr 
rows, Nr columns and characteristic function l(bi^ G C\). 
Since b^^ is a deterministic function of bg, p (bi^lb^) can be 
factorized as follows (FigO: 

p{hR\hs) = p{hR,hQ\hs) 

= p{hR\hQ,hs)p{hQ\hs) 
= l{hReC\)p{hQ\hs) 

The LDGM mapping can either compress or expand the Kr 
quantized bits i.e. the LDGM coding rate can be greater than 1. 
The hR nodes always have degree 2 and they simply perform 
forwarding of messages under the sum-product algorithm. 

E. Scalar Quantizer 

In general, a vector quantizer can be used at R. However, it 
is shown [IJ that QMF performs within bounded gap of capac- 
ity even with a scalar quantizer. Under scalar quantization, the 
observation for every bit from S is quantized independently. 
If each yR^i is quantized into 6Q[Ai] for z = 1,2,... fNs, 
then the p (bg Ib^) function node factorizes into fNs separate 
nodes each representing a scalar quantization operation. 

fNs 

p{hQ\hs) = ]\p{hQ[Ai]\hs,i) , where 

fNs 

U A, = {1,2,...,K^}, A,nA, =0Vz7^j 
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where Ai denotes the subset of indices in bg that observation 
yn^i is quantized into. As a result, the variable nodes of the two 
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Fig. 4. Factor graph: LDPC code at S, LDGM code at R with 1 bit scalar 
quantizer. 

Tanner graphs are connected by function nodes representing 
the stochastic relations p {bQ[Ai]\bs,i) among them. Hence- 
forth, these are called quantize (Q) nodes, as they are induced 
by the quantization procedure at the relay. An example factor 
graph showing the LDPC-LDGM construction is illustrated in 
Fig H] where each symbol observation yR^i is quantized into 
one bit (i.e. Kr = fNs and A[i] = {i}). As shown, there are 
four kinds of nodes in the resulting factor graph: observation 
(OBS) nodes, variable (VAR) nodes, check (CHK) nodes, and 
quantize (Q) nodes. Some VAR nodes in the C| subgraph 
share OBS nodes with VAR nodes in the subgraph. This 
is because of multiple access at D. 

E DBLAST Scheme 



ZSR 




Fig. 5. Channel model using DBLAST 

The factor graph shown in Fig. [4] can be simplified fur- 
ther using the Diagonal Bell Labs Space-Time architecture 
(DBLAST) ||29ll . As discussed in the Appendix, the degree 
2 OBS nodes (representing multiple access) are factorized 
using DBLAST. Under DBLAST, the destination observes two 
orthogonal sets of observations (see Fig. [5]). The factorization 
is shown in equation ([T]) where Y := [Ysd Yrd]- 

An example of the simplified factor graph is depicted in 
Fig. [6l In this graph, VAR nodes in the two Tanner graphs 
are connected only through Q nodes. Since yRD,i = for 
z = 1, . . . , fNs, we rename jRD^ifNs+j) = Yrdj, for j = 

The resulting graph has a structure similar to an irregular 
LDPC code but with special Q constraints. In Sec III-Gl 
sum-product updates for this graph are derived following the 



general principle outlined in 1251 . It is shown that for a simple 
one-bit quantizer each Q node further factorizes into a CHK 
constraint and a dummy VAR node. This reduces the factor 
graph to a Tanner graph that does not have any special nodes. 
Such a property is useful to leverage existing techniques used 
for the design of low-power, high-throughput LDPC decoders. 



LDPC Check Constraints 



LDGM Check Constraints 




i ioBs i i 

f{yi\bs,i) f{y2\bs,2) f{ys,Ns\bs,Ns) 



i i i 

f{yR,l\bR,l) f{yR,Nn\bR,Nn) 



Fig. 6. Simplified factor graph with one-bit scalar quantizer and DBLAST. 

G. Decoding Algorithm 

For the point-to-point system, belief-propagation is an it- 
erative algorithm that computes the a posteriori probability 
to decode message bits. The algorithm computes this exactly 
if the factor graph has no cycles. Otherwise, it computes the 
approximate a posteriori probability for each bit ll25ll . For the 
factor graph in Fig [51 messages being passed on the edges 
of the factor graph and the update rules at the variable/check 
nodes stay unchanged. The only new ingredient in the mix is 
the Q nodes introduced by our framework. 

Let the subscripts V, C, and Q denote VAR nodes, CHK 
nodes, and Q nodes respectively. For F G {C, Q}, let ujyl 
denote the message sent from variable node V to function node 
F in the iteration. Every edge in the graph is connected to 
exactly one variable node and a message on the edge represents 
the a posteriori probability for the respective variable. The 
messages can be represented as LLRs, but for the sake of 
simplicity we represent the messages as a two-dimensional 
vector in this subsectiorQ. cjvf •= [Po Pi], where ujyf{l) = 
Po G [0, 1] represents the probability that the bit is and 
ujyf{2) = pi e [0, 1] represents the probability that the bit is 

1 (Po+Pi = 1). 

The message sent from V to F G {C, Q} is the normalized 
product of all incoming messages into V except for the 
message from F. The normalization ensures that po -\- pi = 1 
for the outgoing message. The message sent from C to V is the 
indicator function that the check is satisfied, marginalized on 
the bit represented by V^ The message sent from Q to V is the 
marginalization of the function p {bQ[Ai]\bs,i) on the symbol 
represented by V. bq is computed from a noisy observation 
of bs,i, the node Q imposes a probabilistic constraint on the 
variables. Since the quantization is scalar: Vu G {0, Ijl^^l and 

^e{0,i}, 

g{u,v) :=p{bQ[Ai] = u\bs,i = v) . 

This function is fully represented by a lookup table with 
2 1 1+1 values, which is used to derive the update rule for 
Q. 



^ Later we will .^j^. 
in belief propagation 



replace uj by w, the commonly used message log ^ (LLR) 
tion. 
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fNs Nr 

/(X|b5,bQ) = Yl f{ySD,i\^S4)Y[f{ySD,ifNs+j)^yRD,ifNs+j)\bQJ^bs,ifNs+j)) 
i=l j=l 
fNs Nr 

= n fiySD,i\bs,i)Y[f{ySDXfNs+j)\bs,ifNs+j))fiyRD,ifNs+j)\^Qj) 
i=l j=l 
Ns Nr 

= Y[f{ySD,i\bs,i)Y[f{yRD,ifNs+j)\bQj)^ 

i=l j=l 



(1) 



[1 - pf pf] jp^ 

\\ LDGM Check Constraints 




f{yi\bs,l) f{y2\bs,2) f{ys,j,^\bs,Ns) fiyRAM fiyR,Nn\bR,Nj 

Fig. 7. Equivalent factor graph of that in Fig.|6] The Q nodes are factorized 
into a CHK node and a dummy variable. 

As an example, let us consider a one-bit scalar quantizer at 
the relay and derive the update rule. For this case, the Q node 
can be further factorized into a CHK node and a dummy VAR 
node that sends a constant message. Note that Ai = {i}^ i = 
1,2,..., fNs and Kr = fNs. The factor graph is depicted 
in Figured Consider Wu G {0, 1} and v G {0, 1}, 

g{u,v) :=p{bQ^i u\bs,i v) 

=(1 — pf)l{u = v} -\-pfl{u 7^ v} 

here Pf \ erfc \J~^^^ denotes the probability of bit error 
for scalar one-bit quantization over a BIAWGN channel. Since 
the function g is symmetric in u and it can be assumed that 
the VAR node is of the source, and the marginalization is on 
V. Let the other VAR node be V . This leads to the following 
update rule: 

^Qv(l) = (1 -P/)cjv'q(1) +p/a;v'Q(2) 
cjQv(2) = (1 -p/)cjv'q(2) +p/cjv'q(1), 

This takes the same form of a CHK node update with incoming 
messages ujyQ and [1—pfPf]. Therefore, the Q node in this 
set-up specializes to a CHK node with additional dummy VAR 
nodes sending constant message [1 — pf pf] that depends on 
SNR^i^. The resulting factor graph is depicted in Fig. [71 

HI. Code Design 

In this section, design of specific codes for QMF relaying 
is discussed. Typically sparse graph codes like LDPC and 
LDGM are drawn randomly from ensembles, which are de- 
scribed using degree profiles. In the point-to-point case, if the 
block-length is sufficiently large, the decoding performance 
of such codes converges to the ensemble average |22|. Let us 
consider degree profiles (A^, ps) and (A^^, pr) for the LDPC 



and LDGM codes at source and relay respectively. A^ and 
ps are polynomials representing variable and check degree 
distributions for the LDPC code: 

oo oo 

^s{x) = ^A5,iX*"\ ps{x) = ^ps,ix'~^ 

i=2 i=2 

Here Xs,i and ps^ denote the fraction of edges with degree i 
at a variable node and at a check node respectively. For the 
LDGM code, we have the similar definition for Ai^ and pR 
except that these are regarding the edges connecting check 
nodes and variable nodes for bg (not b^^). 

These profiles must satisfy the following constraints: 

^^-^ Jo ps{x)dx Kr ^ jl pR{x)dx ^ f 

Xs{x)dx' XR{x)dx 1-/' 

Xs{1)=M^)=Ps{1) = Pr{1) = 1 

In addition to the sub-graphs representing the two component 
codes, the joint factor graph shown in Fig. [7] also includes 
edges connecting them (via Q nodes). As discussed previously 
in Sec III-CI let us consider a fixed one bit scalar quantizer. 
The edges connecting with Q nodes are considered fixed in 
the rest of this section. In contrast, the edges in LDPC and 
LDGM subgraphs are drawn randomly from the ensemble 
using construction procedure described in f2J\ . 

In point-to-point channels, the typical method to analyze and 
design sparse graph codes is to compute the ensemble average 
performance for given degree profiles assuming infinite block 
length (convergence to computation trees). The ensemble aver- 
age performance (decoding error probability for given SNRs) 
is calculated using density evolution developed in ED, ES- 
The two key elements of classical density evolution, namely, 
concentration around ensemble average and convergence to 
computation tree channels for sufficiently large block length 
hold for the proposed QMF relaying system as well. The 
proofs can be readily extended from those of point-to-point 
channels ||22| EB- 

Without loss of generality it is assumed that the all-zero 
codeword is transmitted from S. This is a result of the 
symmetry of the relay channel. However, R does not transmit 
the all-zero codeword because the source-to-relay channel is 
noisy. For sufficiently large block lengths and a given value 
of SNR^i^ there is a typical sequence bg that is mapped 
to a typical hR based on the LDGM code. The probability 
of occurrence for atypical codewords vanishes as the block 
length becomes large. Therefore, it is ignored for computing 
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the ensemble average performance. A typical bg comprises of 
Kr{1 —pf) O's and KRPf I's. pf is defined in Sec. III-GI For 
a given degree profile, (Xr^pr), each bit of the typical b^^ 
is i.i.d. Bernoulli (g'), where q is the probability of having an 
odd number of I's in a column of the generator matrix (drawn 
randomly from the LDGM ensemble). 



E 



pRU)/j A 1 - (1 - ^PfY 



To develop density evolution rules for QMF relaying, we 
consider the belief-propagation algorithm with log-likelihood 
ratios as the messages passed among variable nodes and 
various function nodes. Let w^y and w^l denote the message 
sent from the function node F to the variable node V and vice- 
versa, at the l-th iteration. F G {C5',Ci^,Q,0s',0i^} represent 
the LDPC CHK nodes, LDGM CHK nodes, Q nodes, OBS 
nodes at S and R respectively. V G {V5',VQ,Vi^} represent 
the VAR nodes corresponding to b^, bg and b^^ respectively. 

The sum-product update rules in terms of the commonly 
used LLR's are written as follows: 



(0 

VF 



PGAr(V)\{F} 



(0 

F'V 



(2) 



(0 

ov — 



^v, (0,V) = {(05,V5),(0i^,V^)} 



w^:^^^^ = 2tanh"^ 



tanh 
if F = C5,C^ 




FV 



2tanh"^ | (1 - 2pf) tanh (^^vt) ) 



if F = Q 

Here A/'(-) here denote the set of neighboring nodes and wy 
represents the LLR from channel observation. 

Compared to the point-to-point case where there is only 
one kind of variable node (V) and one kind of CHK node 
(C) the update rules can be expressed simply by ^ and (|3]) 
where F = C. Density evolution analysis tracking the density 
of these messages in each iteration. For the point-to-point case 
with degree distribution (A, p), there is only one type of edge 
and the evolution is expressed using a pair of coupled recursive 
equations as follows: 



i-i) 



Here F (•) denotes a transformation on the density as defined 
in II2TII . (g) denotes the convolution operator and P|^j denotes 

the density of message ^|^.}- represents the conditional 
density of the LLR of the point-to-point channel. 

For the QMF relaying case, there are 4 types of edges and 
densities for messages along all of them must be tracked. The 
recursive density updates are derived similarly: 



Function nodes to variable nodes: 



QVq 



Variable nodes to function nodes: 



Pf 



d(0 
VcCs 



Hi-i) 



(1-/)A5. {P^%,y 

A (^civg) 



d(0 

QVq 



d(0 



=Pvs 



f (1) \^(^) 



(3) <q=Ea..(^Sv«)^ 



Here Sr{-) denotes the Dirac delta function at point r G M. 
6 i-pf shows up in the expressions because the Q node 

is equivalent to a CHK node connected to a constant. The 
differences in evolution rules between the QMF relaying 
and point-to-point channel arise due to the probabilistic Q 
constraints in the joint factor graph. 

As in the point-to-point case, Py^ is the conditional density 
of the LLR of the source to destination channel, given that 
an all-zero codeword is sent from S. Py^ is the marginal 
density of the LLR of the relay to destination channel under 
the marginal law that b^^: is i.i.d. Bernoulli(g). 

The density evolution rules derived above are used to 
compute the probability of error in decoding of b^. For 
successive interference cancellation using DBLAST, b^^ must 
also be reliably decoded. Density evolution rules to compute 
probability of decoding error for b^^ can be similarly derived. 

IV. Bit Interleaved Coded Modulation 

So far, the discussion has focussed on the BMS Gaussian 
relay defined in Sec. III-Al For the high SNR regime, input 
alphabet G A^^ and x^^ G A^^ where A represents 
constellation points in a high-order modulation scheme must 
be considered. In practice many systems use BICM |[2Q|| to 
combine channel codes designed for binary alphabet with 
high-order signal constellations. BICM has also been proposed 
for various cooperative channel scenarios I30II |[3T1 1321 [1161 . In 
this subsection, a procedure is discussed for extending the 
coding framework from the BMS relay channel to a relay 
channel with inputs from high-order alphabets. 
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Under classical BICM 12011 , a point to point Gaussian 
channel is decomposed into parallel independent memoryless 
"sub-channels". Every "sub-channel" PY\B,s{y\^-, ^) has bi- 
nary inputs h G {0, 1} and depends on state s G {1, 2, . . . , L} 
which is chosen uniformly and known to both the terminals 
(2^ is the cardinality of the chosen signal constellation). At the 
receiver, LLR for a bit that was mapped to state s is calculated 
from symbol observation G C (in case of MIMO receiver 



LLR{y, s) = log 



B\Y,S 



{b = 0\y,s) 



PB\Y,s{b = l\y,s) 



However, this binary channel is not guaranteed to be output- 
symmetric i.e. the crossover probability for a bit is not 
independent of its value. Let /a (A) represent the PDF of 
LLR{y,s). The channel is output symmetric if the following 
condition holds: 

/a|b(A|6 = 0) = /a|b(-A|6 = 1) 

Conventional methods for designing linear coding schemes 
such as density evolution etc. cannot be used with asymmetric 
channels. This issue is resolved by adding random dithers at 
every bit to make the channel output- symmetric as proposed 
in |I50|, |[33l, f34l. Dithers are i.i.d. Bernoulli (^) variables 
known to both the transmitter and receiver. For a dither d G 
{0,1} the channel PY\B^s,D{y\b^ s^d) is binary, memoryless 
and symmetric (BMS). 

LLR{y,s,d) = {-l)'^LLR{y, s) 

This method is called parallel BICM (PBICM) in f3?| and 
Fig. [8] shows the architecture for a PBICM point to point link 
having L states i.e. signal constellation of size 2^. {m^}^^ 
represent messages and and b'^ the transmit codewords 
before and after dithering. The equivalent BMS channel can 
be characterized by L, the SNR of the underlying AWGN 
channel and the symbol mapping in modulation. In the rest of 
this paper we consider that Gray mapping is used. 

In order to use PBICM with the relay channel, a definition 
of quantize-and-map operation under PBICM is required. With 
a PBICM modulator at source S, the observations at relay 
R (Yr) represent L interleaved codewords. If R performs 
quantization at the symbol level, then the decomposition into 
independent binary sub-channels will be lost. As an alternative, 
it is proposed that R perform quantization at the bit level. 

S and R both use PBICM modulator blocks with con- 
stellation size 2^ having state and dither vectors given by 
ss^sr^'Ds and D^^ respectively. The QMF operation at R 
is described below (depicted in Fig. O: 

1) For observed symbol sequence ysR '= {ysR,j}j=i per- 
form PBICM demodulation. The output is represented 
as {nsR4}i=i where each n^i^,^ := {nsR,ij}j^i 
represents LLRs for the codeword. 

2) Quantize every LLR in {nsR^i}iLi. As an example, for 
a one bit scalar quantizer this simply involves observing 
the sign of LLRs. 

3) Encode the quantizer output {mR^i}f^-^ using an LDGM 
code. 



4) Transmit the resultant codewords {b^^ i}^^ using a 
PBICM modulator. 
Using this definition of QMF, the Gaussian relay channel is 
decomposed into parallel BMS relay channels. The BMS relay 
channel is shown in Fig.fTO] It is characterized by constellation 
size at S and R and the SNR of the underlying AWGN links 
i.e. SNR5i?,SNR5D,SNR,^z). 
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Fig. 8. PBICM architecture, {d^jf^^ are dithers, and d'- = I - 2di. 
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Fig. 9. QMF relaying with PBICM. 



V. Link Design Example 

In Sec. Unl an extension of density evolution tools was 
developed lEIllEZ I for joint LDPC-LDGM factor graphs based 
on QMF relaying. In this section, a link design example 
with construction of explicit codes is shown for a DBLAST- 
equivalent channel shown in Fig. [5] BMS relay channel. The 
performance of designed codes is presented using simulations 
with high order modulation based on PBICM principles de- 
scribed in Sec. [IVl 

A. System Parameters 

The capacity advantage of cooperative relaying is most 
pronounced when the source to relay link is significantly 
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better than the direct Hnk between source and destination. We 
therefore consider an example scenario where the S to R Hnk 
is 10 dB stronger than the others. 

SNR^D =SNR^z), SNR5i? = 10xSNR5D (4) 

1) Modulation Order: As a guidehne for system design 
use the following information-theoretic bound on maximal 
achievable rate using QMF relaying with continuous Gaussian 
inputs xs and xr and a vector Gaussian quantizer at the noise 
level. 



'^QMF,G = 

r (i-/)(2:g(snr5d) + /(2:g(^ 2 



(5) 



<SNRs 



mm 



SHRsd) 
-f 



Here (Cq (^) •= log {1 -\- x) is the AWGN point-to-point 
capacity at signal-to-noise ratio x. If the inputs are constrained 
to structured constellations such as 16 QAM, 64 QAM, then 
the achievable rate with 2^"^ -QAM modulation and BICM is 
computed as follows: 



'^QMF,n = 

mm < ^ 



(6) 



(1 



/)£„ (SNRsd) + /£„ + SNRsD 

f)€n{SmRD) + €n{SmsD)-f 



Here too we use a vector Gaussian quantizer at the noise level. 
Note that n e {2,3,4} and (tn{x) denotes the 2^'^-QAM 
constellation-constrained point-to-point capacity at signal-to- 
noise ratio X under BICM. 

2) Listening-time Fraction: For QMF, the listening-time 
fraction f at R can be independently optimized to maximize 
system throughput ll35l . lEl, El- The optimal /* is found by 
balancing the two terms in the minimization of (0): 

(1 - n^G (smsD) + r^G ^^^^^^ 



SNR 



SD 



= (1 - n(tG {smRD) + (to (snr^d) - /* 

Alternatively a sub-optimal listening fraction / can be used 
based on reduced channel knowledge at relay. It is shown in 
(2J that this does not have a significant impact on throughput. 

For system parameters in Eq (O, 7^qmf,g and 7^QMF,n 
are plotted for n = 2,3,4 in Fig [TT] vs. SNR^i;). For each 
point, the optimized listening fraction /* is used. To design 
a link with throughput of 5.4 information bits per symbol. 



both 64 QAM and 256 QAM are potentially good choices for 
modulation having QMF information theoretic thresholds at 
14.18 dB and 13.47 dB respectively. Let us choose 64 QAM 
(6 coded bits per symbol) for the example design, which means 
that S should use an LDPC code of rate 7^ = ^ = 0.9. The 
optimal listening fraction corresponding to SNR^i;) = 14.18 



dB is /* 



. This determines the LDGM coding rate 
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Fig. 11. Maximum achievable rate for QMF relaying with modulation 
constraints on channel inputs plotted vs SNR^/^) for SNR relationships in 
Eq. EJ 



B. Code Design 

Codes C| and Cj^ optimized for the above system parame- 
ters can be designed using density evolution tools II2TII . This 
involves finding good degree profiles that have the lowest 
possible decoding SNR threshold and randomly generating 
finite block length codes from them. 

In order to reduce the computational complexity of density 
evolution we use the Gaussian approximation to density evo- 
lution developed in |[36ll . Additionally, we use the following 
heuristics to reduce the search space for profiles. 
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1) For we consider check degree profiles that are 
concentrated |36| i.e. all check degrees (from edge 
perspective) are either /c or + 1 for some integer k >2. 

2) For Cg we consider variable degree profiles with maxi- 
mum degree of 8. 

3) ForC|^ we limit ourselves to regular LDGM profiles. 
Using these heuristics we design the following degree 

profile for the system parameters in this example. 



Xs{x) = 0.2Sx + 0.32x2 + 0.28x^ + O.Ux^ + 0.0009x^ 
ps{x) = OMx 



Xr{x) 



0.96x 

= x' 



Simulation results for the bit error rate in decoding of bo- 
using codes (with block lengths ^ 10^ and ^ 10^) drawn 
from above profiles are shown in Fig. |12(a)| using PBICM 
with 64QAM modulation, one bit scalar quantizer and an ideal 
interleaves As shown the BER performance is < IdB of the 
QMF threshold. For the single relay scenario, the information- 
theoretic thresholds for QMF and CF are identical, therefore 
as a reference for comparison thresholds for DF, AF, and 
the no-cooperation case are also shown. The DF and the 
AF thresholds are computed using the following expressions. 
Derivations follow standard analysis of the schemes and are 
omitted here. 

7^DF,n = max 

/e[o,i] 

mm{f(tn {Smsn) , (1 - f)(tn (SNR^z,) + (t^ (SNR^d)} 



7^AF,n = -^n (SNR^d) 



1 



SNR 



<tn (SNReff) 

SNR5i?SNRi^i) 



SD 



SNRe„ ......^ , i + sNR5K + SNR^^ 

The optimal listening time for DF is determined by the channel 
parameters, while that for AF is always 1/2. 

For the DBLAST architecture, b^^ must also be reliably de- 
coded at or below the target SNR (for successive interference 
cancellation to work). Fig. |12(b)| shows the BER for b^^ which 
is also within < IdB of the QMF threshold for both of the 
block-lengths. 

VI. Conclusions 

The QMF relaying scheme has the following key advantages 
over other known relaying schemes such as AF, DF, and CF. 

1) For the single relay network, it outperforms AF and DF 
at high SNR. 

2) For the single relay network, it achieves the same per- 
formance as CF but reduces channel feedback overhead. 
Unlike CF, QMF does not require knowledge of forward 
channel strength at the relay. 

3) For arbitrary relay networks with multiple relays, QMF 
achieves better high SNR performance than AF, DF and 
CR 

In this paper, a low-complexity channel coding framework is 
developed for QMF relaying. For the single relay network, the 
framework performs within (0.5 — l)dB of fundamental limits. 

The techniques presented here can be extended to complex 
system scenarios, which are discussed below. 
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Fig. 12. Code design simulation results. 



A. Multiple Relays 

When there is more than one relay in the system, the 
proposed factor graph extends in a straightforward manner. 
Optimal listening schedules can be computed for each of 
the relays. As proposed, the source would use an LDPC 
code and each relay would use an LDGM code based on 
its respective schedule. The joint factor graph would include 
multiple LDGM sub-graphs. 

The DBLAST architecture proposed in this paper extends 
naturally to networks with one level of multiple non-interfering 
relays e.g. the diamond network. As discussed previously, 
DBLAST significantly reduces the complexity of the factor 
graph. DBLAST requires that all codewords from relays are 
decoded correctly at destination in order to permit successive 
interference cancellation. This additional constraint does not 
lead to a reduction in the QMF information-theoretic achiev- 
able rate. In fact, such a requirement is explicitly considered 
in the probability of error analysis for the QMF scheme in 0). 
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However, some challenges for multiple relay networks re- 
main to be addressed. When the relays can hear one another 
or the source can reach the destination via multiple hops, 
it is unclear how the DBLAST architecture can be applied. 
In such scenarios, an alternate space-time architecture must 
be considered. Moreover as the number of relays increase, 
the channel knowledge overhead required to compute optimal 
listening schedules becomes large. Practical techniques at 
the physical and MAC layers are required to address this 
complexity. These are considered as directions for future work. 

B. Rate Adaptation and Hybrid ARQ 

In the link design example, suitable coding rates, constel- 
lation and listening fraction are computed for a given set 
of operating channel conditions. However, optimizing codes 
based on instantaneous channel conditions is not feasible in 
practice. Under commonly used rate adaptation mechanisms, 
terminals switch between a few candidate codes and a few 
candidate constellations based on channel conditions. Coop- 
erative links need to consider multiple channel parameters to 
determine transmission rates i.e. for a single relay three SNR 
parameters are required as opposed to just one for a point- 
to-point link. This makes rate adaptation schedules for relay 
networks more complex. An advantage of QMF relaying is 
that rate adaptation schedules depend only on the ability of 
the destination to decode as opposed to DF, where adaptation 
must consider decoding at relays as well. 

Modern adaptation mechanisms like hybrid automatic re- 
peat request (HARQ) can be incorporated into the proposed 
framework. Additional parity bits for refinement sent from the 
source after receiving a repeat request from the destination. 
It can be cooperatively delivered to the destination using 
QMF relaying. The joint decoding factor graph is expanded 
to incorporate these refinement parity bits and the decoding 
algorithm remains unchanged. 
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Appendix 

The QMF relaying scheme introduces correlation between 
xs- and x^^, which can be thought of as coding across transmit 
antennas in a MIMO channel. A natural space-time architec- 
ture for such a channel is DBLAST. Using DBLAST for the 
relay channel has also been proposed in |[T2lll[Tnill29lll3Qll . It 
relies on introducing a delay of one block at the relay and using 
successive interference cancellation (SIC) at the destination. 
At the k-th block the destination receives the superposition of 
the following: 

• signal from the source containing the codeword sent at 
block k, namely, X5'(m/c) 



• signal from the relay containing the side information 
about the source's codeword at block k — 1, namely, 

Xi?((7/c-i) 

Messages sent from the source are independent across blocks. 
At the k-th block, the destination jointly decodes block k — 1 
(message rrik-i and side information Xi^(g/e_i)) by treating 
^^(m/e) as Gaussian noise. The receiver subtracts relay's 
codeword :siR{qk-i) from its received signal Y[k] and keeps 
the residual Y [k] for decoding the next block. This architecture 
allows the use of a simplified equivalent channel model. Note 
that the one-block delay introduced at R has the added benefit 
of allowing time for QMF processing at R. 

1) Simplified Channel Model: The equivalent channel 
model is shown in Fig.[5l For decoding the block ^fc—l message 
m/e_i, the decoder takes two inputs Y[k] and Y[k — 1]. We 
can think of Y[k] and Y[/c — 1] as two orthogonal links with 
independent Gaussian noise. Therefore, for the purpose of 
code design we can alternatively investigate a simpler model 
depicted in Figure [51 In this model, 

Y,, = hij^i + Zij, {ij) = {R,D), {S,D), 

ySR = ^5i?X5 + ZSR 

As an example, let us consider a scenario where D has 
two receive antennas (M = 2). In that case, the DBLAST 
equivalent channel becomes IITtII : 

Yij = hij:>Ci ^Zij, {ij) = {R,D),{S,D), 
where. 



hsD = llhill, huD = ^\\h2±l\\^ + ^^p^\\i^^\\2 

h2±i and h2||i denote the perpendicular and parallel com- 
ponents of h2 with respect to hi, respectively. The signal- 
to-noise ratios of the three links are SNR^i^ = {hsnl'^Ps, 
SMRsD = \hsD\'^Ps, and SNRrd = \hRD?PR respectively. 

Remark 1: Consider the original channel and the DBLAST- 
equivalent channel. Note that the capacities of these two 
channels are within two bits of each other. This is based on 
the following observations: 

1) The min-cut upper bound for both channels are within 
one bit of each other (for any listening fraction f G 

[o,i]j. 

The mutual information across cut {5'}, {R^ D} remains 
unchanged between the two channels. Consider the 
mutual information across the cut {S^R}j{D}. It is 
known that SIC achieves the sum capacity of multiple- 
access channels. In the original channel (Fig. [2]) S and 
R have unlimited cooperation. As a result, the min-cut 
bound for DBLAST incurs a power-gain loss of at most 
(1 - /) bits. 

2) QMF relaying scheme achieves the min-cut upper bound 
to within one bit for the two channels yj. 
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